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ON INTERMEDIATE ORBITS. 


By Dr. G. W. Hrxix, West Nyack, N. Y. 


The assumption of the Keplerian ellipse as the first approximation to the 
motion of a planet leads to some inconveniences, the worst of which, perhaps, 
is that the mean longitude at the epoch suffers a perturbation proportional to 
the time, and thus the mean motion in longitude is not the same in the perturbed 
and unperturbed orbits. In the case of Saturn this perturbation amounts to 
110” per year. This inconvenience would be avoided if, developing the per- 
turbative function in a series proceeding according to the cosines of multiples 
of the angle // contained at the sun by the radii of the planets thus 


R= gy (7,7) + (7, 7’) cos H+ ¢, (7, cos 2// + ..., 


we should suppose ¢, (7, 7’) annexed to the potential of the attraction exerted 
by the sun and the resulting differential equations integrated and the expres- 
sions of the coordinates thus deduced regarded as the first approximation. 
Then the following approximations could be obtained by the method of varia- 
tion of the elements or otherwise by attributing to 7? the value 


R= ¢,(r, 7’) cos // + ¢,(7, cos 2// +.... 


In this way of approaching the question, the first approximation to the 
motion of a system of planets revolving about their central body would involve 
a potential function containing the radii of the planets alone, without their 
longitudes or latitudes ; thus the force acting on each planet would be directed 
along its radius, but would be a function of all the radii. The orbits are then 
all plane curves whose planes pass through the central body, and thus we know 
precisely what functions the latitudes and the reduced longitudes are of the 
orbit longitudes. Also there is equable description of areas by the radii, and 
when the latter are known functions of the time the determination of the orbit 
longitudes is reduced to quadratures. Thus the question is narrowed to the 
finding of the radii as functions of the time. 


Let 7, 7’, ... v, v, ... denote the radii and orbit longitudes of a system of 
planets, m, m’, ... their masses and .V the mass of the central body ; 7 = +m, 
=M-+m,...; the perturbative functions for m,m’,.... Then 
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2 HILL. ON INTERMEDIATE ORBITS. 


the differential equations for determining the radii are 


d*p ‘ dy 
dt’ dt’ or’ 
d*f 
dt dt or’ 


But, 4, 4’, ... being constants, we have 


dv h div 
dt dt 7?’ 


The functions 72, 72’, ... are to be reduced to 
R = n'¢(r,r’) + m’¢(r, 7”) 4 


r) + me 7’) 


where we have the relation ¢ (7, y) == ¢(y,#). To evaluate the function ¢, let 
the symbol m denote the operation of taking the arithmetico-geometrical mean, 


then 
1 1 


¢ (7, 7") M (7 + 7, — M Vr? — 2)’ 


where it is supposed that 7 > 7. The two quantities under the symbol ® may 
be brought still nearer to equality by continuing the operation, and writing 


=a, yr* —r =f, we have 
1 1 
And if 7 and 7’ are not too near each other a sufficiently approximate value is 
2; 2 


(7, 7’) - 4 ° 
(a + 3) af(a 4+ 
If we write 


1 


2 r 2 xr? 


the differential equations determining 7, 7”, . .. become 


d*r oY oY 
Vit > — 
dt or’ dt or? 
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HILL. ON INTERMEDIATE ORBITS. 


From which we derive the integral equation 


dr? 
dt dt ( 


C’ being the arbitrary constant. The left member of this cannot be a negative 
quantity. Consequently, if we construct in a space of » dimensions, x being 
the number of planets, the surface whose equation is 


04+c=0, 


the representative point 7, whose coordinates are the values of the variables 
r,7,..., must lie on the positive side of this surface. In such a system of 
planets as our solar system composed of the eight major planets, at least one 
fold of this surface is closed ; and it is within this that the representative point 
? always falls. Thus, in this case, we are able to set definite inferior and 
superior limits to 7, 7’, ..., which, however, in general are not the min/ma- 
minimorum or maxima-maximorum values of the variables, since the point /’ 
only attains the surface when all the planets are together on the lines of their 
apsides. 

A particular solution of the foregoing system of differential equations can 
be obtained in the following way: Put 


20 20 
and solve these equations, regarding 7, 7’, ... as the unknowns. Let 
, 
ra, 


be a system of values satisfying them. Then the last equations evidently sat- 
isfy the differential equations. When this process is applied to the solar 
system the values found for the a@ do not greatly differ from the mean distances 
of the several planets from the sun. The excessive smallness of the mm rela- 
tively to V/ brings it about that there is but one real solution ; and quite ap- 
proximately, by neglecting in 2 all the terms dependent on the functions ¢, we 
get 


a= 


By substituting these values in the terms arising from these functions we can 
arrive at more exact values for the a; and this process can be repeated as 
often as is deemed necessary. Evidently all this comes to the adjustment of 
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1 HILL. ON INTERMEDIATE ORBITS. 


the planets at such distances from the sun that for the given values of the m 
and the / the centripetal force may be equal to the centrifugal force in each 
case, thus rendering it possible that all the orbits may be circles. 

Desiring now to find inferior and superior limits to the particular variable 
7", in the surface 2 — ( = 0 we suppose this variable to be a function of all 
the rest, and thus that when it arrives at a maximum or minimum value we have 


except that ~ , = 1. In order to find limits for the values of 7” we take the 
group of equations 
oY 
= 
20 
and in it replace the equation .~ = 0 by the equation 2+ (=0. Solving 


these equations, regarding 7, 7’, ... as the unknowns, we shall always, in cases 
like our solar system, be able to find two sets of corresponding real values for 
these quantities such that the two values of 7 shall contain between them the 
observed value of this radius from which the value of C was derived. And no 
other roots can be found which satisfy this condition. In ordinary cases these 
roots are non-multiple and there is but one solution to the problem. 

For the sake of illustration, limiting ourselves to two planets, we may take 
a hypothetical case, suggested by Jupiter and Saturn. Assume as the values 
of the planetary masses 


1047.355 ’ 


Wt = 


~ 3501.6" 
From Leverrier’s Tables for the time 1875, Jan. 1.0, we get 
log 7 = 0.7367630 , log = 0.9953504 ; 


and, with the Julian year as the unit of time, 


log = 8.0841634n, log — 8.876073. 


dt dt 


From the same source the values of the osculating elements necessary for the 
computations are 


log a = 0.7162505, loge = 8.6865018, log a’ = 0.9794334, loge’ = 8.7291462. 
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HILL. ON INTERMEDIATE ORBITS. 


Accommodated to the units of length and time we employ, 
logu = 1.5967563, log = 1.5964658, log = 8.5762479, log m’ = 8.0520753. 
log A? = 2.3119809, log h” == 2.5746499. 


When these numerical values are substituted in the equation 2 + C' = 0, it is 


found that C = — 0.166513246. Employing brackets to denote numbers cor- 
responding to common logarithms, the equation 2 4+ (= 0 becomes in this 
case 
(0.1730042] — [0.5871988] + [9.6485411] — [0.3256952] 
+- [6.6283282] |” | = 0.166513246 . 
us 
, oY 

And the equations 0, = 0 become 

2. [0.5871988] — [0.1730042] + [6.6283232] 0, 


2 [0.3256952] — [9.6485411] — (6.6283232] E |= 0, 
| J 


where ¢ denotes the same function as Laplace’s 4" is of a. 

The solution of the first and third of these equations gives inferior and 
superior limits for 7, the solution of the first and second the same for 7’. Em- 
ploying the tentative process we obtain as the corresponding values of the 


logarithms of the radii in the four solutions 
log 0.6926819, 0.7390339, 0.7152497, 0.7152334; 


0.9776183 0.9238410, 1.0390527. 


log r — O0.97T76575 


The solution of the second and third equations gives the values of the 
radii for which the planets could describe uniformly circles. These are 
log a = 0.7152396 , log a = 0.9776403. 
If we assume that the maximum eccentricities of Jupiter and Saturn are 


respectively 0.0594902 and 0.0845677, the minimum and maximum values of 
log 7 and log 7’ are severally 


log r = 0.6895907, log r = 0.7413343, log 7’ = 0.9411223, log r’ — 1.0147527. 


The differences between the former values and these are, of course, due to 
the neglect of the terms of /2 involving the angle /7/. 
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For the purpose of graphically exhibiting the numerical results just ob- 
tained we may suppose a system of rectangular axes to be drawn with 7 and 
r as coordinates. The four sets of corresponiling values of these variables 
give four points on the oval which constitutes one of the branches of the curve 
having 2 + C= 0 as its equation. These suffice to roughly indicate the course 
of the oval within which the representative point /’, exhibiting the simulta- 
neous values of the radii, always keeps. The position for which both radii are 
constant is well towards the centre of the oval, If we draw the right line 
bisecting the right angle at the origin formed by the axes of coordinates, the 
intersection of this line or its nonintersection with the oval will show whether 
there is intrusion or not of the spheres of the planets on each other. In the 
case here treated there is no intrusion ; each planet maintains its character as 
exterior or interior ; consequently the distance between them can never vanish, 
and stability of motion is assured. 

The given example shows that there must exist a large domain in which it 
is possible to expand the function #2 in a series of powers and products of the 
deviations of the radii from their mean values. Here we will limit the exposi- 
tion to the case of two planets, as the formule given can be readily extended 
when there are more. Let us put 


where the second terms are small relatively to the first. The differential equa- 


tions then become 
dr oY =o 
dt? ar dt? 


In developing 2 in a series of ascending powers and products of x and « we 
will stop with terms of three dimensions. ‘The constant term can be omitted 
as its retention would only serve to modify the value of the arbitrary constant 


20 20 
As aand are corresponding roots of the equations ~~ 0, = @, the 
cr 


terms of one dimension vanish. Also the term in #2’ can be supposed absent, 
for, if present, it could be removed by a linear and orthogonal transformation 
of variables. We write therefore 


2 aa? — + ca? + + Var? + 4 


where «, a’, b, b’, c, & are constants.* 


* The reader is asked not to confound these a@and @ with the quantities previously denoted by 


the same letters. 
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HILL. ON INTERMEDIATE ORBITS. 7 


Adopting two constants ~ and ~’ at present left indeterminate, but to be 
determined hereafter so as to fulfil certain conditions, the differential equations 
can be written 


+ = — a?) + + + 
dt 

AP 


Before integrating these it will be well to ascertain superior and inferior 
limits for z and #’. Those for x are determined by the equations 


C aa? — + ca* + bate’ + + 4 


— +- Cr vax’ L hy? 


If we put, = + ! “C’ this will be a small quantity of the same order as x 
a 
and the limiting values of « can be developed in an infinite series of powers 
of 7. Thus 


c= 94 34 7 aa 1 Gat gg 61% 
3 2 aa” 18 a aa ata” 27 a 


where we take 7 positively for the superior limit and negatively for the inferior. 
The corresponding value of z’ can be derived from the equation 


a a’ 
he limits for 2’, by putting 7 = = *~, can be obtained from those of . by 
a 


removing the accent from the symbols which there have it and applying it to 
those which are destitute of it. 

By neglecting the right-hand members of the differential equations as last 
written their integrals are simply 


2=ecosy, =—é€cos¢’, 


-and é being arbitrary constants and ¢ and ¢’ are arguments increasing pro- 
dy’ 


By substituting these expressions 
( 


portionably to so that 
( 


for 7 and a’, 
cx? + Qhra’ = (ce? + Ve?) + 4 ce’ cos 2¢ + bee’ cos(¢ + ¢’) 


+ hee’ cos(¢g — ¢’) + Ve" cos 2¢’, 
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8 HILL. ON INTERMEDIATE ORBITS. 
cx + +- ba? = 4 (ce? + be?) + 4 ce? cos 2y’ + cos (¢ + ¢’) 
+ Bee’ cos (y — ¢’) + be’ cos 2¢. 


Whence it follows that z contains the additional terms of the second order 
with respect to ¢ and ¢’ 


1 9 1” 1 2 hee’ 
—, (ce + — ce’ cos COS(Y + ¢’ 
Fn? | —(n + +¢) 
hee 1 
n?— (n — (¢ —¢) 2 n* — 4n” 


and x’ terms which are obtained from those of « by simply interchanging the 
accents. 

Pushing the approximation to terms of the third order, 2? is found to con- 
tain the following terms of this order: 


a? = [5 ce? + cos ¢ + COS — —, cos 3¢ 
nr 4n? — n® 6 n° 


+ COS (2¢ p’) 
n= — (n + Cr+ + n° — (n — 


" On? — 4n? cos + + 3 — cos (p — 2¢’). 


The similar terms of z* are obtained from those of «* by interchanging the 
accents. The terms of the third order of xz’ are 


1 + be’ h(e*—1é 1 ce + We” 
zr’ = | at cos + 9 s— + ( cos 


2 4 — n2 


1 be’ 1 
cos 3¢ 


b’ 

: - cos (2¢ + ¢’) + ; —, — cos (2¢ — ¢’ 

E +n én | 2n on’ —n 6ni Bn (4¢ 


, cos (¢ + 2g’) + ,cos (y — 2y’). 
E Ln ] 2 (¢ ¢) on — 2n (¢ 


cos 


n 6n’ 
The coefficients n> — a and n* — a? are of the second order, hence limiting 


ourselves to terms of the third order 


(n? — = (n?—a@)ecosg, (n* —a”)a 


= — a”) cos ¢’. 
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The constants e and ¢’ being arbitrary, the coefficients of cos ¢ in v and of 
cos ¢’ in # need receive no corrections in the following approximations. Hence 
the coefficient of cos ¢ in the right member of the first differential equation 
should vanish, as also that of cos ¢’ in the second. This gives us the two 
equations 


of which the second is obtained from the first by interchange of accents. 'To 
the degree of approximation to which they are pushed they serve to determine 
the values of » and as functions of the constants a, a’, 4, ¢, As 
the next terms which would appear in these equations through the following 
approximations are of the fourth order with reference to end ¢’, the errors of 
n and nv’, as determined from these equations, are of the same order. 

The remaining terms of the third order in the right member of the first 
differential equation are 


+ 1 — cos 3¢ +14 |e” cos 3¢ 
* — 4n° n — 
; ee cos(2¢ + ¢’) 
be , be bb’ bh’ 
+ ; » | cos(2¢ — ¢’) 
— n™ as — — n* 
4n™ — n* Zan’ n° 
Ve he’ b” 2 
4n™ — n* — n’ 


The similar terms of the right member of the second differential equation are 
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obtained by interchange of accents in the preceding expression. The terms 
of the third order in « are then 


2nn Qnn’' * 4n? — 1 3n? + 4nn 
he be hh’ hb’ ee (2 
F Ae’ /? ee” L Dep’ 
| ; » Cos (¢ + 


The similar terms in are obtained from these by interchange of accents. 
On many occasions it will be interesting to know the equivalent of C, the 


arbitrary constant attached to the integral equation, in terms of «, a’, 4,4’, ¢, &, 


e,e. To obtain this we compute the non-periodic terms of 2”, ye eve 
dt 
The non-periodic terms of 2°, yer can be derived from these by inter- 
Af 


change of accents. Limiting ourselves to these terms we have 


(4n° 4n- 4n” 


Substituting these values in the integral equation and eliminating 7° and 7”? by 
means of the two equations which determine them, we find 


l6a™ (4a* — 144 a° 
Si! + 5a",. , 8a*— + 5a‘, 
2u “(40° a“) (4a a’y a* 


16a* (4a? — a*)* 144 a” 
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It is necessary to notice the effect of the vanishing of the divisors intro- 
duced by the integration. In the terms of the second order if ~ — 2” the 
coefficient of cos 2¢’ in « and the coefficient of cos (¢ — ¢’) in «’ apparently 
become infinite. Also, if ”’ = 2v, the coefficient of cos (¢ ¢’) in w and the 
coefficient of cos 2¢ in «’ are in like case. This evidently is the analytical 
warning that in these cases the arguments 2¢, 2¢’, ¢ — ¢’ cannot be distin- 
guished from ¢ and ¢’. Thus the coefficients of their cosines should have 
some indetermination. But let us see what conditions the quantities a, a’, b, 0’, 
c, ¢, e, © must satisfy that either of these relations may have place. We shall 
not allow that any of the quantities a, a’, 7,’ can be infinite. Then, from the 
equations determining « and ~’, it is plain that, in the first case, either 4’ — 0, 
or — ¢* = 0, and in the second case, either — 0, or & = 0. Thus 
all the mentioned coefficients in # and 7’, instead of becoming infinite in these 
cases, really take the form °. 

The particular supposition of ¢ — 0, e’ — 0 makes the first approximation 
vanish, and it is plain from the differential equations that 0, x 0 form 
a particular solution of them in which and »’ are quite indeterminate. This 
‘ase therefore does not demand further consideration. But suppose that, in 
the first case, 4’ = 0 or, in the second, / — 0. Then, in 2, the term 4x2” or 
the term /7*2’ disappears, and the terms in z and # having indeterminate co- 
efficients should also disappear. However, the relation ~ 2n’ or n’ = 2n, 
without being exactly fulfilled, may be so very nearly. Let us take the first. 
If, in the equations determining ” and we make = and = except 
in the first terms and those which involve the divisor 4" »#*, we shall have 
quadratic equations for the determination of “° and «*. Thus, putting 


a? — a? — (8 c? + 35) © — 4 be’ + 2 B) 

= + (fc? + b'o— phe’ — 

(4 Ae 13 c*) 

the solution of these quadratic equations will give 


n® = 4a” + 4) — 6*(2e* — e”), 


4n’- n= — 28 2) (2e- 


The ambiguity of the sign before the radical is determined from the following 
considerations. When 7 is a somewhat large positive quantity it is evident 
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that we ought to take the upper sign, as always approximately ”* = «*, and 
then n° = 40° + 4,3, thus #* is larger than 4°. Hence, if we begin with 7° 
larger than 47° and keep diminishing ,%, continuity demands that the upper 
sign be maintained. And we may persist in diminishing / until, if 2e° — e* is 
a positive quantity, +-) 6°(2e — e*) when n*= a? — 41 6° (26 — e*) and 
n° = 407; thus 77 is still larger than 47°. But if ,4 is a somewhat large nega- 
tive quantity it is evident that the lower sign must be attributed to the radical ; 
and then approximately #* = @*, n* = 4a° + 43; that is, here 7° is less than 
4n°. Then may be numerically diminished until 6° (2e — e*), when 
(28 — &*) and and is still less than 42°. On 
the other hand if 2¢* — ¢* is a negative quantity, 4 in both cases can be nu- 
merically diminished until it vanishes, and then, in the first case, we have 


4) (é* Ze"), n® = +- 2¢"), and, in the second 
case, = — hy (e* — Qe*), = — 4y (e* — 2e*). 

From all this it is apparent that we cannot have 41% — 1° = 0 unless two 
conditions are satisfied. The first is 3 0, and the second either 4’ 0 or 
2 -*- 0. Tf the second condition is not fulfilled the limits towards which 
tend the coefticients = and severally belonging to cos 

Zn 4n™ (n n'y 
2¢ in w and cos(¢ — ¢) in #’, can be written so as to cover all the cases 
and 


The ambiguous sign within the radical must be taken so as to render the quan- 
tity following positive ; that without the radical depends at once on the sign 
of 4 and on whether ,t has been supposed to approach the limit from positive 
or negative values. These limits are independent of 4, and it will be perceived 
that instead of being of the second order with respect to e and é they are of 
the same order as the first terms of « and a’. The second supposition of the 
second condition viz. 2c? — ¢* = 0 renders these limits infinite ; but it is prob- 
able that this conclusion would be modified if terms of the fourth order were 
included in the equations determining ~ and w’. 

There is no need of giving the similar investigation for the assumption 
n — 2n. All the results of the latter are obtained from those of the former 
assumption by simply interchanging the accents. Thus the limits of the co- 
efficient of cos 2¢ in.’ and the coefficient of cos (¢ — ¢’) in « are severally 


4, + (2e? — + (2e? — 
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In all this we have assumed that the quantity under the radical sign in 
the expressions for 7° and ~” must not be negative. This is simply to insure 
stability of motion. If the constants of the problem are such as to make it 
negative, or if either »* and x” come out negative, the problem has still a real 
solution but the motion is unstable. 

The question of periodic solutions is intimately connected with the van- 


ishing of integrating divisors. If we suppose that /n ~ vx’ — 0, where / and 
? are integers prime to each other, by adopting an argument ¢/ such that 
d¢ dy d¢ dé : , . 

— 7 and , the series for 2 and 2 reduce to the form 

dt dt at dt 


vor a, cos 4+ 3, sind + a, cos 3, sin 


The variables then return to the same values after ¢/ has gone through a cir- 
cumference as also do their differentia! coefticients. If we conceive a system of 
rectangular coordinates for graphically exhibiting the values of » and ’, the 
representative point 7’, in its motion, will deseribe a certain curve which, if the 
variables are confined within limits, will lie within a limited portion of the 
plane. This curve may have multiple points ; let us see what is the condition 
necessary and sufficient that this may have place. The motion of # and 2’ 
being a constant oscillation between minimum and maximum values, the curve 
must be either a spiral having no multiple points, in which case there would 
be a progressive diminution or augmentation of the minimum and maximum 
values of 2 and «’, or there must be an infinite number of multiple points. The 
periodic solution is the case which lies between these two classes of solutions. 
We may conceive the angle of intersection at the multiple point to constantly 
diminish through variation of the arbitrary constants of the problem until it 
vanishes, when the periodic solution is reached. 
Let us suppose that the solution of the system of differential equations 

dt on’ df 
admitting the integral 


has been discovered and that it may be written 


Let ¢ denote another value of the independent variable ¢; then, if the equa- 
tions 
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¢and ¢ being regarded as the unknowns, are satistied by real values of ¢ and 7¢, 
the curve under discussion will have multiple points ; but, if the roots of these 
equations are all imaginary, the curve will be a spiral. It appears that when 
v and a are contained within finite limits the roots ¢ and ¢ are all real and 
infinite in number ; and thus that, in general, there are an infinite number of 
multiple points. 

Let us now suppose that the arbitrary constants are adjusted so that, in 
addition to the equations f, (¢) = f,(¢). £,(¢) = f, (2) being satisfied by certain 
values of ¢ and /, the equation 

(4) 


is satisfied by the same values. (Accents have been used to denote differen- 
tiation of the form of f.) Then it is evident that the last equation is the con- 
dition of the presence of a periodic solution ; and if we put ¢ — ¢=— /, the 


equations 


fit+P)=£(), £(¢+ P)=£(0, 


are satisfied for all values of ¢. 

A periodic solution must have at least one arbitrary constant less than in 
the general case. But, for the differential equations we are treating, there is 
usually a periodic solution in which there are only two arbitrary constants. 
This happens when both variables arrive simultaneously at their maximum or 
minimum values. Integrating the differential equations by the application of 
Maclaurin’s Theorem, we should have, in general, 


it ht eee, 

D 


ce 


where a, a, 4, 4, are the arbitrary constants, and in the derivatives of 2 it is 
understood that and are to be replaced severally by @ and a. In the case 
where both variables arrive simultaneously at a maximum or minimum value, 
and we count ¢ from this epoch, not only 4 = 0, 4’ = 0, but also all the coefficients 
of the odd powers of ¢ in « and 2’, as is easy to see from the successive differ- 
entiation of the derivatives of 2 with respect to the time. We may therefore 
put 


a’,t* + +.... f, (¢). 
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In the previous equations for determining the times of the representative point 
/ passing through multiple points, we may suppose ¢ ft, and the equations 
f,(— t) = f, (0), f,(— 2) = f,(¢) are fulfilled for all values of 7. Also the equa- 
tion showing the presence of a periodic solution is, in this case, an identity. 
This is explained by the circumstance that the point /’? having attained the 
limiting curve whose equation is 2 — ( = 0 returns on its previous path. 
Hence the condition necessary for the existence of a periodic solution in this 
vase is that when the path is sufficiently prolonged backward it should again 
meet the limiting curve for real velocities. The motion of /’ will then be a 
swinging back and forth on this are. This is tantamount to saying that the 
ralues of « and a must be so adjusted that a value ¢, for ¢ can be found which 
will make « and « have such values that they satisfy the equation 2 + (= 0. 

It may be of interest to know the earlier terms of the expansions of « and 
«in powers of #. For the sake of brevity we put 


(3a)? 
[ 
Then 
x ( x ‘ 
r l or l 2 . 3 . 4 
on | E cr | 46! 
== 


1.2 ow’ 1.2.3.4 


( dai Oren ox”? on on 6! 
where in the partial derivatives it is supposed that 2 and «’ are replaced sev- 
erally by @ and 

Evidently, for these series proceeding according to ascending powers of 
#, we may substitute others proceeding according to cosines of multiples of an 
angle proportional to ¢, or 7 and ” may have the form 


vor a =a, + 4, cos nt - a, cos 2nt 


But here we are not entitled to assume that the coefficients rapidly diminish. 
To bring about such a relation as x = n’ it may be necessary for the arbitrary 
constants we have denoted by e and e’ to be of the zero order of magnitude. 
If in the expansion of 2 according to powers and products of . and « we 
omit all the terms involving both variables the integration of the differential 
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equations is reduced to quadratures. Let us write in a more general way this 
expansion of 2 
It must be noticed that the coefficients of the terms involving both ~ and a’ 
are of the order of the planetary masses while those which multiply the simple 
powers of or « are of the order of the central mass. If we put 
2 
&, yy C yy 
+ + aya 4 
the differential equations resulting from this expression of the potential fune- 
tion are 
dt 
From these are immediately derived the integrals 
da 
da? 
and thence the final integrals 
di 
JV (CH + aye + +...’ 
( dx 
c= . 
JyV(C + + +...) 


The polynomials under the radical sign, on being equated to 0, will have each 
two real roots, which, on account of the smallness of « and 2’, will be, sever- 


ally, in the neighborhood of the values 
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The exact values of these can readily be found by trial ; let them be denoted 
severally, by a, b, a, b’. Then we can write 


da 
A 
Jy {la r) b) (A,+ A gt” ) 
dx 


The limiting values of 7 are a and b, and those of 2’, a’ and b. Within these 
limits the last factors under the radical sign are always positive and do not 
vary much, and thus x and ~@ are periodic functions of ¢, the periods being 
given by the formule 


yp 4+ Ay’ +...) 
Li’ 
] [(a x’) b ) (A 0 A T A 
b 


By adopting, in place of x and «’, variables ¢° and ¢ such that 


= (a b) (a b) cos 22’ = + b’) (a’ b’) cos  , 


the integrations are considerably simplified. The motion of ¢/ and ¢ is unlim- 
ited like that of 4 Let 


1 
[A,+ B, cosy 4+ 2B, cos 


[A’, + 4+ 4+... = B cosy’ + 22, cos 2¢ 


Then 
sing + B,sin®’ +..., 


+ Bsin + B,sin 2’ 


The values of ¢/ and ” can be derived from these transcendental equations 
and have the forms 


g=n(t+c)+ FE, cosn(t+c) + F,cos2(nt+c)4+..., 


¢)+ cosn' (t+ ¢) 4+ Ey cos 2 (n't + 4 
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From these we can pass to the values of cos ¢ and cos ¢”, and thence to those 


of w and «’, which have the form 


t= 4 + COS ¢ COS 2¢ y 
=4G,+G/ cos¢’ + G, cos 2y’ + ..., 
‘ 


y and ¢’ being put for ¢ and vt 
The reversion of the periodic series can be avoided by the use of definite 
integrals, and we can pass directly to the values of the G, and G/. For we 


have 


us 


with a similar formula for G;. By adopting ¢° as the independent variable 


this becomes 


») 
[(a + b) — (a — b) cos ¢| cos [4 7B, sin ¢ + 7B, sin | 


In the particular case of 7 0, this reduces to 


G a+b b) B,. 


Proposing now to take into account the omitted terms of 2, which we will 
denote as 2, let the values of # and .’, just determined, be denoted as x, and 
vw, their corrections as dr and dv’, the differential equations for the latter will 
be 


Pox 
») 1 Sn | 


These can be integrated by successive approximations, introducing into the 
right members the values of # and 2’ from the previous approximation, begin- 
ning with #-~ 7, 2 = a’, In order to prevent ¢ from getting outside of the 


functional sign cos, we adopt the following device: Let ” and x’ be two in- 
determinate constants, whose values come out only at the end of the approxi- 
mation, but which are equivalent to the rigorous rates of motion of the 
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arguments ¢ and ¢’. We add n*dx to each member of the first equation and 
n*dx to each member of the second. Then, putting 


= [n? + + Say (e+ + (2? + + 2,7) +...) de 


and Y for the similar quantity in the second equation, we have 


dt dt 


In any stage of the approximation we shall have 


Y= cos(¢e +e’), cos(i¢ 4 ve’), 
from which 
d2 = — (in + cos (ig “¢’), 
= (n+ cos (ig + vy’). 


Then the equations for determining « and w’ are 


The terms of 2 in *”’ and xx? can be removed from it by a linear trans- 
formation, but at the expense of reintroducing a term in vz’, as also a term 
involving the product of the velocities into the expression for the living force. 
For sabstituting 

z=y ly, y ly, 


if /and / are so determined that 


3a, + dy (1 + 2d’) + (2 + W) + 3a,l? = 0, 
Baul? + (2 + Ul’) + ay (1 + Ql’) + = 0 
the terms in 2 involving y°y and yy* vanish. And if we put 
== 
B,, = 2ayl + a, (1 + + Qa 


? 2 
B, = + a,l + 4 


a,t + ay”, 


v2 
9 
> 3 2 
we have 
2= By + Byyy + Bay? + + Buy® 
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and the expression of the living force becomes 


a, ay dy dy » dy? 


and the differential equations of the problem are 


» 
(1+ + 0) = + + Buy - 


The elimination of / or / between the equations which determine them 
vives rise to an equation of the fifth degree, which consequently, always has a 
real root. As, and «@, are of the order of the masses of the planets, ap- 
proximate values for / and /’ are / 

In integrating the preceding equations one may follow either of two 
courses ; first, we may neglect the terms in 7’ and y”, which will afterwards 
be regarded as perturbative, and then we shall have linear equations with con- 
stant coefficients easily integrable ; or second, we may regard the last terms of 
both members of each equation as perturbative, and the equations restricted 


to the remaining terms are integrable by quadratures, as has been shown. 
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PROOF OF A FORMULA DUE TO CAUCHY. 
By Pror. W. H. Ecuons, Charlottesville, Va. 


The function 


(a?ar)" 1 a ( 2\n 
vanishes when = (¢ = 0,1,..., p being any finite quantity. 


If we expand (1) with respect to the first column, we have 
1 1y (1 “ nt+r—1) a n) 
(a —1)...(a" —1) 


This rational integral function of « of the mth degree has the » roots a~”, 
a’’,...,a"'?; whence we have the identity 
t=n-—1 a n+r—1) ] 


Evidently, if « > 1, we have 


[=e 1)" Pt! 
— = 1+ 3 an”. 
=1 (a—1)...(a’ 1)" (3) 


r 
Let « = «* ”, where ¢ is any positive integer or zero, then 


=1(a@ —1)...(a" —1) 


In (2) let p — 1, we then have Cauchy’s formula (Comptes Rendus, 
1840. Chrystal’s Algebra, II, 316), and (3) becomes Euler’s theorem (Introd. 
in Anal. Inf., § 306). 

As a consequence the sum of the products, taken ¢ at a time without repe- 
tition, of the quantities 


(a—1)...(a*— 
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ON SOME FORMS OF LAGRANGE’S INTERPOLATION FORMULA. 
By Pror. W. H. Ecuors, Charlottesville, Va. 
1. The well known form 


a, 

’ ( ) 
C7 (a, . 


which has for its value 


R= (x (x7 


(v lying between the greatest and least of w, @,,....«,) may be written out 
according to the first row, in the standard form 


a;) 


Ear (2) 


or by the first row or column according as we give precedence to the + sign of 
vor p, in the formula 


3) 


wherein /?,— 1, /’, is the sum of the products of the quantities 1/ 
at a time, without repetition. 


taken p 


The values of 
Ia, 
(a,— a,)’ 


for the three cases «, r,t, = 1, a, = a’, are, respectively, 


= (1) 


ri(n—r)) 


(ii) 


i 
n 
r= 
Ai 
it 
i 
/ 
(1 7 ] nN). 
p ] ’ 
(1 1/n)...(1 r/n) 
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which becomes, when x. 


Np ’ 


‘%» being Fourier’s number of the 7th order and class p. 


1)... (a 1)(1 L/u)...(1—1 an (111) 
in which 
(— 
q=0 te" qa 1) 
Sp a’ 


8 (a 1)... (a? 1). 


/ 


The corresponding values of /? are in these cases 


R= (— 1)" (1 —a#/a)... (1 — 2/a") at (iii) 


2. Instead of writing the determinant (1) in solid columns, we may write 


it 
0 ’ ’ an (4) 
ni n ’ 
ay 
wherein 
= —a,)...(@#—a,). 
(4) gives 
n 1 no 
fax 2 A, fa, - . (9 
p= we + 1)! ) 


in which we have for computing the A’s, A, = 2"'/w,"', and 
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3. Again, we have 

lo 2) 
wherein 
(6) gives 
n > 
& fru, (7) 
f=} 
a in which we have for computing the A’s, A, == fa,/a,'", and tL 


+ (—1)"' fa,. 


CHARLOTTESVILLE, July, 1893. 
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SOLUTIONS. 


330 


Finp the sum of the series 


P4+3?+ 6410+ 157+...+ [an(n 4+ 1)f. 
[| Artemas Martin.] 
SOLUTION. 


The fifth differences vanish and the first terms of the difference series are 
1, 8, 19, 18, 6, 0; hence 


C,.8 +,C,.19+ ,C,.18 
where ,,C,, ete. are the binomial coefficients. [ Geo. R. Dean.) 
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Four equianharmonic points give four triangles which have four cireum- 
circles. Show that the inverses of any point with regard to these four circles 
are equianharmonic. [Frank Morley. | 

SOLUTION, 


It is proved in Salmon’s Conic Sections, § 54, that if the perpendiculars 
from three points a,, @,, @, on the sides of a triangle ¢,c,c, meet at a point a, 
then also the perpendiculars from ¢,, ¢,, ¢, on the sides of a,7,a, meet at a 
point ¢,. And this can be readily proved by elementary geometry, by observ- 
ing that the condition of either concurrence is, that in the hexagon a,, ¢,, 4, 
C,, @, ¢, the sum of the squares of three alternate sides is equal to the sum of 
the squares of the other three sides. 


Now 
Co — Cs SIN SIN 4,4,4, 
: 
SIN €,C,C, SIN 
¢—¢,|_ sinegc, sina,a,a, 
C, — ¢, sina,a,a,’ 
whence 
Cy— _ Sin d,aa,. sin aat, a, — a,.a, — a, 
C, — C3. SIN . SIN a, — a, 


Hence two corresponding cross-ratios of the two tetrads of points have equal 
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: absolute values. They have also congruent amplitudes ; for 
am(c,—¢,) am(a,—a,)—7/2, 


am (c,—¢,)  am(a, — a,) + 7/2; 
and, therefore, 
am (c, —¢,.c,— ¢,) am(4, —a,.a,—a,). 


Hence corresponding cross-ratios are equal. 

In particular, let ¢,, ¢,, ¢, be the centres of the circles @,4,4,, 
then c, is the centre of the circle a,¢,a,; and the theorem is: The inverses of 
the point x , with regard to the four circles through three of any four points 

: a, have the same cross-ratios as the four points. This is a covariantive state- 
ment, in which we can substitute any other point for the point o. 
| Frank Morley.} 
339 


Snow that the areas of the curves 
y= pe’ +gqe+r, and y= + pr+qr+r, 
taken between the limits 2 + 4 and x — A, are given by the formule 
Q = 2h + Rph®, and 2 = Bh + (Fp + 


respectively. [| W. Echols.) 
SOLUTION I. 


From the well-known formula 


< we have for the first curve 
Q= (pat + ge + r)de pa + + 
h 
2A ( px? + ga + + = + ph’ ; 
“J and for the second curve , 
r+h 
= = + pa? + ga + r)dx =} mat + pe? 
= h 
| = 2h + pe + ge+r)+ 2inx) Qyh + 2mzx) h>. 
[Artemas Martin.]} 
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SOLUTION II. 
By the familiar formula called Simpson’s rule the mean value of 
y= me + 
is the sixth part of the sum of the extreme values plus four times the middle 


value. But 
(e+ AP + (2 — hf + 42° — Go’ + GA*z, 


(2 Ay + (a — + = 62° + 2h’, 
(29 +h) +(x2—h) + 62; 
Yn +- pa + qr-+r 4 ph? 
y+ +h p), 
and A = 2hy + p). 
Putting m — 0, we have the other result. (W. ML. Thornton. | 


SOLUTION IL. 
As applications of 


r-h 


we have for the first one, 
A 3 
- 
2=hy 4 
for the second 


¢ ‘ 
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A WoopEN hemisphere floats in water, vertex down, with 1/n of its axis 
immersed. Find the specific gravity of the hemisphere. 
[Artemas Martin.| 
SOLUTION. 
Let s be the specific gravity of the hemisphere and , its radius. The 
weight of the hemisphere equals the weight of the water displaced, or, 


r ) 
4 
which gives 
Sa — 


[ W. O. Whitescarver.]| 


| 
‘4 
3 
| 
ji 
| 
| 
| 
4] 
| 
| 
} 
y 
| 
3 
= 
Py Jos 
an. 
4 
a 
= 
F 


I8 SOLU'TIONS—EXERCISE. 
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A sotip sphere and a solid cylinder of equal radii roll from rest down the 
same inclined plane ; compare the times of their descent. 
| Artemas Murtin. | 
SOLUTION. 


By the principle of vis viva, 
mv? + — 2mgs sini. 


The motion being pure rolling, 


fe. 
For the sphere, 
[= mr . 
For the cylinder, 


Substituting these values and reducing, we have, for the sphere, 


as 
dt gs Sin 7, 


and for the cylinder, 


= -= gesini. 
dt 
Integrating, we get 

whence 

t 

15° 


[S. 7. Moreland; M. C. Andrews; CU. Mendenhall; R. Dean.) 


EXERCISE. 


357 
GIVEN tan,w = ty, wherein 
tan,w = x — + 


w=-etw, 
and 
x= m(cos 7 sin 
determine wu and v as real quantities in terms of x, y, m, and 3. 
[Jrving Stringham. | 
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